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Abstract 


Regenerative  liquid  propellant  guns  (RLPGs)  have  been  studied  for  many  yeaK.  RIPG 
firings  almost  always  show  large,  high-frequency  pressure 

phenomenon,  a  fluid  dynamics  model  of  the  combustion  chamber/gun  tube  of  an  RLPG  has 
developed.  High-frequency  oscillations  are  generated  naturally  by  die  physics  modeied  m  Ae 
code.  Dimensional  analysis  has  been  applied  to  the  governing  equations  to  furnish  insight  into 
the  dynamic  relationships  among  pressure,  density,  temperature,  droplet  size,  etc.  I  e 
dimensional  analysis  predicts  a  scaling  law  between  different-sized  fixtures.  ,The  co^p^  !^ 
results  exactly  match  these  predictions.  This  agreement  establishes  the  basis  of  self-sinulanty 
in  the  physical  processes  (e.g.,  tradeoff  between  droplet  size  and  burning  rate)  and,  m  addition, 
helps  verify  the  computational  model. 
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1.  Introduction 


Pressure  oscillations  associated  with  combustion  chamber  instabilities  in  regenerative  liquid 
propellant  guns  (RLPGs)  have  been  studied  for  many  years  to  attempt  to  identify  sources  of  the 
instability  and  to  develop  techniques  to  mitigate  the  pressure  excursions  [1,  2].  Over  the  past 
year  and  a  half,  the  U.S.  Army  Research  Laboratory  (ARL)  has  been  exploring  the  use  of 
“dimensional  analysis  theory”  to  simplify  the  description  of  RLPG  combustion  processes  and  to 
provide  a  firm  theoretical  foundation  for  research.  Dimensional  analysis  theory  can  significantly 
reduce  the  number  of  key  variables  that  need  to  be  analyzed  in  order  to  efficiently  search  for 
ways  for  mitigating  combustion  chamber  instabilities.  This  report  summarizes  the  results  of  the 
study. 

A  diagram  of  the  30-mm  concept  VIC  (6-C)  RLPG  is  shown  in  Figure  1.  The 
monopropellant  in  the  liquid  reservoir  is  prepressurized  and  located  between  the  control  piston 
and  the  injection  piston.  An  igniter  injects  hot  gas  into  the  combustion  chamber.  As  the 
chamber  is  pressurized,  the  control  piston  is  pushed  to  the  left,  opening  the  injection  orifice.  The 
motion  of  the  control  piston  depends  on  the  balance  of  forces  acting  on  the  piston  and  includes 
the  forces  resulting  from  the  pressure  generated  in  the  chamber,  the  retarding  forces  acting  on  the 
reservoir  side  of  the  piston,  and  the  retarding  forces  associated  with  the  damper  assembly.  The 
injection  piston  follows  the  control  piston,  injecting  the  propellant  from  the  reservoir  through  the 
annular  orifice.  The  combustion  takes  place  in  the  combustion  chamber,  and  the  gas  then  flows 
into  the  gun  tube.  There  is  a  large  area  change  from  the  chamber  to  the  tube. 

Gun  firings  have  been  simulated  using  a  lumped  parameter  code  RLPGUN  [3-6].  The 
injected  liquid  is  assumed  to  instantaneously  break  up  into  droplets  and  ignite.  The  droplet 
diameter  is  obtained  from  an  experimental  correlation  involving  the  entrance  velocity,  gas 
density,  liquid  density,  etc.  [7].  Although  the  functional  form  appears  relevant,  it  is  necessary  to 
multiply  the  droplet  diameter  obtained  from  this  correlation  by  a  constant,  which  depends  on  the 
specific  injection  system,  to  obtain  the  proper  pressure  rise  rate.  This  constant  may  depend  on 
droplet  coalescence;  liquid  in  the  intact  core;  a  droplet  ignition  delay;  or  other  simplifications  in 
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Figure  1.  ARL  30-nun  Concept  VIC  RLPG. 


the  breakup  model,  which  do  not  fully  represent  the  physical  phenomena.  The  discrepancy  could 
also  be  due  to  extrapolating  data  from  low  pressures  in  the  experimental  correlation  to  very  high 
pressures  with  higher  spray  densities.  However,  using  this  modified  correlation,  good  agreement 
with  data  taken  in  the  tests  has  been  obtained  for  a  wide  variety  of  RLPG  firings. 

Essentially,  all  firings  of  RLPGs  show  large-magnitude,  high-frequency  pressure  oscillations 
that  cannot  be  modeled  using  a  lumped  parameter  code.  To  consider  spatial  variations,  a 
two-dimensional/three-dimensional  (2-D/3-D)  model  LPOSC  has  been  written  for  the 
combustion  chamber/gun  tube  of  an  RLPG  [8,  9].  The  piston  motions,  injection  rate,  and 
injected  droplet  size  are  all  obtained  from  a  lumped  parameter  code  simulation.  The  propellant 
combustion  rate  is  simulated  based  on  an  engineering  correlation.  Due  to  lack  of  information 
about  spray  behavior  under  gun  conditions,  the  very  simple  breakup  model  and  combustion  rate 
equation  from  the  lumped  parameter  code  have  also  been  used  in  LPOSC.  Pressure  waves  are 
generated  naturally  by  the  model,  and  reasonable  comparison  with  experimental  data  has  been 
demonstrated. 

Dimensional  analysis  is  a  mathematical  tool  that  rigorously  allows  the  most  efficient  way  of 
establishing  the  minimal  set  of  dimensionless  functions  that  are  required  to  analyze  and  interpret 
a  physical  problem  such  as  the  RLPG  [10-13].  In  addition,  it  permits  the  construction  of 
experiments  that  are  properly  scaled.  These  scaled  experiments  can  be  used  to  efficiently 
organize  and  interpret  the  results  of  experiments  and  computer  simulations  that  have  already 
been  performed.  Additionally,  they  can  be  used  to  verify  the  internal  consistency  of 
mathematical  models  used  in  computer  simulations  [14]. 

We  have  taken  the  first  several  steps  toward  identifying  the  key  parameter(s)  responsible  for 
combustion  instabilities  using  the  dimensional  analysis  approach.  These  findings  have  provided 
complementary  insight  that  has  been  accrued  by  ARL  and  other  organizations  into  the  nature  of 
combustion  instabilities.  For  example,  pressure  oscillations  do  not  scale  between  the  30-mm  and 
155-mm  guns.  Dimensional  analysis  provides  insight  into  the  physical  reasons  for  the 
differences.  As  a  natural  consequence  of  our  studies,  we  have  also  used  the  scaling  aspect  of 


3 


dimensional  analysis  theory  to  verify  the  mathematical  correctness  of  gas  and  liquid  transport 
developed  in  the  ARL  models. 


2.  Dimensional  Analysis  Considerations 


2.1  Basic  Concepts.  Dimensional  analysis  has  probably  been  used  by  nearly  all  scientists  at 
one  point  or  another  in  conjunction  with  simplifying  a  differential  equation.  Consider  for 
example,  the  differential  equation  for  the  velocity,  v, 

m— +0v3  =  0,  (1) 

dt  F 

with  the  initial  condition,  v(t  =  0)  =  v0.  The  analytical  solution  of  equation  (1)  is,  of  course 
trivial,  but,  if  we  could  not  solve  it  in  closed  form,  we  would  write  the  solution  as 

v  =  v(m,0,vo,t).  (2) 

As  observed,  v  appears  to  be  a  function  of  three  parameters  and  time.  Introducing  the 
dimensionless  velocity  and  time  variables,  y  =  v/v0  and  x  =  fiv^t/m,  respectively,  reduces 
equation  (2)  to 


— +  y3  =  0,  (3) 

dx 

with  initial  condition,  y  =  l.  The  universal  solution  of  equation  (3)  is  y  =  y(x).  A  single 
solution  describes  the  entire  parameter  space. 

Comparison  of  equations  (1)  and  (3)  indicates  the  significant  reduction  in  dimensionality  that 
is  possible  by  making  appropriate  substitutions  to  render  the  equations  dimensionless.  However, 
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is  it  possible  to  reduce  dimensionality  if  we  do  not  know  the  differential  equations  but  do  know 
the  input  parameters  and  variables?  The  answer  is,  “Yes,”  and  it  is  made  possible  by  the  classic 
work  of  Buckingham  [10]  that  was  developed  in  1915.  Buckingham  showed  how  to  reduce  the 
number  of  parameters  and  variables;  he  called  the  reduced  set  “II-functions.” 

The  concepts  of  dimensional  analysis  are  routinely  used  to  improve  our  understanding  of 
complex  problems  in  the  field  of  fluid  mechanics  and  to  organize  the  experimental  data  in  the 
most  efficient  way  (e.g.,  the  use  of  Reynolds  number).  In  some  instances,  the  analytical  form  of 
the  solution  has  been  found  from  purely  dimensional  arguments  (e.g.,  the  blast  wave  equation 
[11]).  A  review  of  these  formal  methods  is  rendered  in  the  next  section. 

We  have  explored  the  use  of  Buckingham’s  theorem  in  selected  subsets  of  the  liquid  rocket 
and  RLPG  problems  (see  sections  3.2  and  3.3).  The  approach  generated  too  many  II-functions 
to  be  of  any  short-term  use.  We  were,  however,  able  to  mathematically  demonstrate  the 
equivalence  between  this  approach  and  the  direct  approach  of  working  with  the  differential 
equations.  Despite  its  lack  of  practical  utility  in  this  study,  the  formal  theory  may  be  useful  in 
future  investigations. 

2.2  Review  of  Dimensional  Analysis  and  Self-Similarity.  The  basis  of  Buckingham’s 
IT-theorem  is  summarized  as  follows.  Assume  a  system  can  be  described  by  “n”  quantities, 
which  can  either  be  physical  constants  (e.g.,  universal  gas  constant),  independent  variables  (e.g., 
time),  or  dependent  variables  (e.g.,  velocity,  temperature),  and  “m”  dimensions.  The  dimensions 
used  in  fluid  mechanics  that  are  relevant  to  our  study  are  mass  (M),  length  (L),  time  (T),  and 
temperature  (0). 

Now,  let  Uj,u2,u3...un  denote  the  aforementioned  n  quantities.  The  solution  of  the  problem 
yields  a  functional  relationship  of  the  form 

F(Uj,u2,u3...un)  =  0.  (4) 
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If  ni,n2,n3...etc.,  are  dimensionless  groupings  of  the  u’s,  Buckingham’s  theorem  proves  that  a 
relationship  of  the  form 

f(n1,n2,n3...nn_m)=o  (5) 

can  be  derived.  The  dimensionality  of  the  problem  has,  in  principle,  been  reduced  from  n  to 
n  -  m.  The  variables  of  the  problem  are  the  n  -  m  II-functions.  The  reduction  from  n  to  n  -  m 
variables  can  be  very  significant  if  n  is  not  too  large.  For  example,  if  n  =  6  and  m  =  4,  the 
problem  is  reduced  to  only  two  variables.  In  such  cases,  the  solution  may  be  accomplished 
analytically  or  by  the  solution  of  two  differential  equations.  On  the  other  hand,  if  n  =  15  and 
m  =  4,  the  number  of  variables  is  1 1  and  not  much  is  gained. 

Equation  (5)  represents  n  -  m  equations  in  n  unknowns.  The  solution  of  such  a  system  is 
achieved  using  established  matrix  techniques  of  linear  algebra  and  leads  to  n  -  m  relationships 
between  the  n  unknowns.  For  tutorial,  as  well  as  practical  purposes,  we  dispense  with  the  matrix 
formalism  and  present  the  methodology  based  on  the  straightforward  procedure  developed  by 
Streeter  [13]  for  determining  the  ITs. 

•  Step  1:  Select  m  of  the  u’s,  which  collectively  contain  m  of  the  dimensions.  For  example, 
if  M,  L,  and  T  are  the  dimensions,  we  must  select  three  of  the  u’s  that  collectively  contain 
M,  L,  and  T.  Therefore  m  =  3  in  this  case. 

•  Step  2:  Label  the  u’s  in  the  manner  consistent  with  step  1.  Using  m  =  3  as  an  example,  we 
select  upu2,  and  u3to  be  the  three  quantities  that  collectively  contain  the  dimensions  M,  L, 

and  T.  We  then  label  u4,u5...un  as  the  remaining  variables  or  parameters. 

•  Step  3:  The  n  -  3  II-functions  are  then  constructed  as  follows: 
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IIj  =Ui‘U^!U31U4 

n2=u?ujpu?us 

nn-3  =u^n'3u2“'3u3n'3un-  (6) 

•  Step  4:  The  final  step  in  the  method  is  to  express  the  11’s  in  their  dimensions  and  then 
determine  the  x’s,  y’s,  and  z’s  so  that  all  the  ITs  are  dimensionless. 

The  aforementioned  procedure  is  most  readily  explained  by  the  example  of  a  point  explosion 
in  air.  This  is  the  classic  problem  of  a  point  explosion  that  has  formed  the  foundation  of  many 
basic  concepts  of  shock  wave  theory.  An  extensive  treatise  on  the  subject  from  the  dimensional 
viewpoint  can  be  found  in  Sedov’s  book  [11].  It  is  presented  here  to  show  the  recipe  for 
constructing  the  dimensionless  II-functions. 

For  the  point  explosion  in  air,  there  are  four  variables:  the  radius  of  the  shock  wave,  R;  the 
density  of  the  ambient  air,  p;  the  time  after  the  explosion,  t;  and  the  energy  of  the  detonation,  E. 
The  dimensions  of  these  quantities  are 

.  R:  L, 

•  p:  ML"3, 


•  t:  T,  and 


.  E:  ML 2T2. 
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The  selection  of  variables  is 


U]  =  R,  u2  =  t,  u3  =  p,  and  u4  =  E.  (7) 

Because  there  are  only  four  variables,  there  is  only  one  II-function.  It  is  determined  by 
inserting  equation  (7)  into  equation  (6)  using  the  dimensionality  of  the  variables.  We  have 

Ili  =  RxtypzE  =  Lx Ty (Mz  L'3z)  (MI^T2) .  (8) 

Setting  the  exponents  of  L,  T,  and  M  equal  to  0  gives  the  following  coupled  set  of  equations: 

x-3z  +  2  =  0,  y-2  =  0,  and  z  +  1  =  0 .  (9) 

The  solution  is  y  =  2,  z  =  - 1,  and  x  =  -  5.  Substituting  the  result  into  equation  (8)  then  gives  the 
dimensionless  function 


IIj  =  R"5t2p"lE .  (10) 

In  section  3.3,  we  construct  sets  of  II-functions  for  the  lumped  parameter  model  of  the  RLPG 
using  the  foregoing  technique. 

2.3  Differential  Equation  Approach  to  Scaling  of  Burning  Rate.  In  section  2.1,  we  have 
shown  how  to  construct  a  dimensionless  differential  equation  from  a  dimensional  equation  by 
recognizing  and  using  the  inherent  dimensionality  of  obvious  physical  quantities.  This  method  is 
not  only  useful  but  is  frequently  the  preferred  method  of  constructing  dimensionless  equations 
that  exhibit  universal  behavior.  The  major  caveat  in  using  this  approach  is  that  our  mathematical 
model  must  be  correct.  If  the  physical  description  is  not  correctly  described  by  the  equations, 
erroneous  scaling  relations  and  predictions  are  possible. 
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The  formal  theory  of  dimensional  analysis  does  not  require  a  knowledge  of  the  systems’ 
differential  or  integral  equations  but  only  that  we  know  what  variables  and  parameters  are 
involved.  On  the  other  hand,  it  provides  insight  into  the  relationship  between  the  variables  and 
parameters.  Combined  with  a  physical  model  of  a  system,  dimensional  analysis  can  provide 
guidance  and  proper  scaling  for  experiments  and  simulations. 

We  call  the  method  of  constructing  dimensionless  equations  from  the  original  equations  the 
“direct  assault  approach.”  The  purpose  of  the  analytical  exercise  of  this  section  is  to  show  how 
this  method  is  used  to  explore  the  types  of  dimensionless  variables  that  arise  with  the  liquid 
propellant  burning-rate  expression.  The  reduction  in  the  dimensionality  of  the  problem  and  the 
ensuing  scaling  relationship  are  developed. 

In  this  section,  the  combustion  process  is  examined.  A  simplified  model  problem  is 
considered.  In  this  model,  the  volume,  V,  is  constant.  There  are  initially  a  specified  number  of 
droplets,  all  with  the  same  diameter.  The  total  mass  in  the  system  (gas  plus  liquid)  is  constant. 
The  dynamic  process  is  the  conversion  of  the  liquid  into  gas,  releasing  energy  and  increasing  the 
pressure  in  the  volume.  The  number  of  droplets  remains  constant  as  the  droplets  shrink  in  size. 

Using  the  assumption  that  the  total  mass  is  constant,  there  is  a  relationship  between  the  gas 
density  and  the  droplet  size.  The  bum  rate  is  a  function  of  the  pressure.  The  process  begins  at 
t  =  0 .  The  initial  conditions  are  the  volume  of  the  system,  the  diameter  and  number  of  droplets, 
and  the  density  and  temperature  of  the  gas.  The  initial  pressure  is  obtained  from  the  Noble- Abel 
equation  of  state.  In  mathematical  terms,  the  droplets  are  allowed  to  reach  zero  mass.  This 
allows  the  total  number  of  droplets  to  be  conserved  in  the  theoretical  model. 

The  density  is  the  average  gas  density  per  unit  volume,  that  is,  the  gas  mass  divided  by  the 
volume.  However,  we  must  account  for  the  fact  that,  when  there  are  droplets  within  the  volume, 
the  actual  gas  density,  po,  will  be  much  higher.  We  now  calculate  this  increase.  The  21 
“primary”  parameters  and  variables  are  defined  in  Table  1.  What  we  mean  by  primary  is  that 
any  other  parameters  and  variables  can  be  expressed  in  terms  of  this  primary  set. 
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Table  1.  Primary  Set  of  Parameters  and  Variables 


Parameters/ V  ariables 

Definition 

V 

Volume 

b 

Covolume 

Po 

Initial  Gas  Pressure 

po  j 

Initial  Average  Gas  Density 

p 

Instantaneous  Gas  Pressure 

p 

Instantaneous  Average  Gas  Density 

N 

Number  of  Droplets  Contained  in  Volume,  V  (Remains  Constant) 

fDO 

Initial  Radius  of  droplets 

Td 

Instantaneous  Radius  of  Droplets 

[  R 

Universal  Gas  Constant 

TG 

Gas  Temperature 

Tgo 

Initial  Gas  Temperature 

Chemical  Energy  Per  Unit  Mass  in  Liquid 

A 

Burning-Rate  Pre-Exponential 

B 

Burning-Rate  Exponential 

Cm 

Specific  Heat  Per  Unit  Mass  for  Gas  at  Constant  Volume 

Cp  =  Cm  +  R 

Specific  Heat  Per  Unit  Mass  for  Gas  at  Constant  Pressure 

Cy  —  PCjrCnn 

Specific  Heat  Per  Unit  Volume 

c 

Ratio  of  Specific  Heats 

Pl 

Intrinsic  Density  of  Liquid 

t 

Time 

Using  the  foregoing  definitions,  we  have 


(ID 


and 


10 


(12) 


R 

7-1* 


The  internal  energy  per  unit  mass  for  the  gas  is 


u  =  CmTG) 


(13) 


and  the  enthalpy  per  unit  mass  for  the  gas  is 


h  =  u  H - —  u  +  RTg  +  pb  =  (Cm  +  R)Tg  +  pb  =  C  +  pb . 

Pg 


(14) 


There  are  two  equivalent  forms  for  the  equation  of  state: 


_  _  PGRTG 
l-bpG 


(15a) 


and 


Pg  ~ 


P 

RTg  +  pb 


(15b) 


We  must  account  for  the  fact  that,  when  there  are  droplets  within  the  volume,  the  actual  gas 
density,  pG,  will  be  much  higher.  In  Table  2,  we  show  the  relationship  between  p  and  pG  with 
the  help  of  additional  parameters  and  variables  that  are  derived  from  the  primary  set. 
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Table  2.  Additional  Parameters  and  Variables  Derived  From  the  Primary  Set 


Parameters/ V  ariables 

Definition 

N 

n  =  — 

V 

Number  Density  of  Droplets  (Remains  Constant) 

4k  3 
vD=TrD 

Instantaneous  Volume  of  Droplet 

o 

i? 

II 

o 

Q 

Initial  Volume  of  Droplet 

'HdO  ”  *^D0 

Initial  Fraction  of  Volume  Occupied  by  the  Liquid 

1— ^Ido 

Initial  Fraction  of  Volume  Occupied  by  the  Gas 

T1d  =  ™d 

Instantaneous  Fraction  of  Volume  Occupied  by  the  Liquid 

Instantaneous  Fraction  of  Volume  Occupied  by  the  Gas 

Pgo  "  1  Pn 

1  ^Ido 

Initial  Tme  Gas  Density 

p.-/ 

Instantaneous  True  Gas  Density 

1  T1d 

mD0  =Pl'udo 

Initial  Mass  of  Droplet 

mD  =Pl^D 

Instantaneous  Mass  of  Droplet 

Mlo  =  NmD0 

Initial  Total  Mass  of  Liquid 

Ml  =  NmD 

Instantaneous  Total  Mass  of  Liquid 

MG0  =  Pgo(^  _  ^Ido)^ 

Initial  Mass  of  Gas 

Mg=Pg(1-tid)V 

Instantaneous  Mass  of  Gas 

m  =  mlo+mgo=ml+mg 

Total  Mass  in  System  =  Constant 

The  first  step  in  the  analysis  is  to  establish  the  relationship  between  the  dimensionless  droplet 
size  and  the  dimensionless  true  gas  density.  The  starting  point  is  the  mass  conservation  equation. 


M  =  Mlo  +  Mgo=Ml+Mg  (16) 


and 


NmD0  +pGo(l-ilDo)V  =  NmD  +pG(l-riD)V 


(17) 
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Dividing  by  the  volume,  V,  and  making  other  substitutions  gives 


npLvD0  +PGo(1-nVDo)  =i1Plvd  +pG(l-m)D). 


We  now  introduce  the  following  dimensionless  variables 


and 


go  = 


Pgo 

Pl  ’ 


Equation  (18)  at  first  becomes 


r|D0  (1  -  x)  +  g0  (1  -  riD0 )  =  yg0  (1  -  xt)do  ) , 


which  then  gives 


T)DO(l-x)  +  g0(l--nDO) 

§o(l  —  x11do  ) 


Here,  r|D0  and  g0  are  constants  that  depend  only  on  the  initial  conditions. 


y  =  y(x,go’'nDo)- 


(18) 

(19) 

(20) 

(21) 

(22) 

(23) 

For  brevity,  we  write 

(24) 
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The  burning-rate  equation  is 


rD  =ido  -  AjpB(t')dt'. 

0 


(25) 


Defining  a  dimensionless  radius  by  the  equation 


Jt-3- 


(26) 


1D0 


and  inserting  the  result  into  equation  (25)  gives 


X  =  l-<t»(t), 


(27) 


where 


4>(t)  =  —  ipB(t')df. 

rD0  0 


(28) 


We  then  have 


x  = 


_Hd_ 

VD0 


(29) 


and 


pB(t)  = 


XD0 

A  dt 


d<|)  _ 


rDo  dA, 
A  dt  ' 


(30) 
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Using  equation  (15)  applied  for  the  initial  conditions  gives  the  following  form  of  the 
dimensionless  equation  of  state: 


P  _  y(x,go,TlDo)(l-bpGO) 

(  rp  ^ 

lG 

_  y^go^ooXl-ao) 

f  rV  ^ 

ag 

Po  l-y(x,g0,r|D0)bpG0 

Joo  J 

i-y(x,g0^Do)a0 

T 

lG0  J 

(31) 


where 


a0=bpG0.  (32) 

In  order  to  define  a  dimensionless  time  and  determine  the  dimensionless  time  dependence  of 
the  pressure  and  other  quantities,  it  is  necessary  to  express  (TG/TG0)  in  terms  of  the 
dimensionless  droplet  radius.  The  required  expression  is  now  developed. 

When  volume  of  droplet  changes  from  UD0  to  uD ,  the  total  amount  of  energy  liberated  by  all 
the  droplets  is 


E  =  eLpLN(i)D0  -  vD )  =  eLpLNvD0  (1  -  x) . 
where  eL  is  the  propellant  chemical  energy  per  unit  mass. 

The  temperature  in  the  gas  is  computed  from  the  equation 


ATr,  —  Tr_  Tgo  — 


eLpLN^D0(l-X) 


CM, 


CmMG 


where  E  is  given  by  equation  (33)  and  the  gas  mass  is  given  by 


Mg=Pg(1-tid)V. 


(33) 


(34) 


(35) 
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Using  the  equation 


where 

K  Mdo(1-x)  (40) 

T  ('nDo(i-x)+go(i-'nDo)) 

The  first  step  in  determining  a  time-dependent  equation  for  the  pressure  is  obtained  by 
inserting  equation  (40)  into  equation  (31).  We  have 
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(41) 


V_ 

Po 


y(x,go,r|DO)(l-a0) 

l-y(x,g0,TiDO)a0 


(l  +  Kt(Po’^1do’§o4  x)). 


Equation  (41)  shows  that  the  time  dependence  of  the  pressure  is  determined  by  the  time 
dependence  of  the  dimensionless  radius,  X  =  x1/3 ;  that  is,  p(t)  is  an  implicit  function  of  time 
through  the  relationship p(t)  =p(A,(t)).  The  other  factors  in  equation  (41),  g0,P0>TlDo>  and a0 , 
are  constants. 


The  explicit  time  dependence  of  the  pressure  is  determined  from  equation  (30).  We  write 


dk 

dt 


ApB 


App  yB(?t3,g0,r|DO)(l-a0)B 
rDo  (l-y(A.3,g0,riD0)a0f 


(l  +  K-r^Ti^go,!-!3))8. 


(42) 


Introducing  the  dimensionless  time 


T  = 


(43) 


where 


T  = 

■*“D0 


1D0 

Apo 


(44) 


then  gives  the  total  dimensionless  equation 


dk 

di 


=  -F(X,Q0), 


(45) 


where 


(46) 


F(A,,Q0) 


yB(A.3>go>'nDoXi-ao)B 

(l-y(X,3,g0,r|D0)aof 


(l  +  KT(P0,TiD0,g0,l -X3)J 


and  Q0  is  the  set  of  dimensionless  parameters:  a0,r|D0,g0,  and  P0 . 


The  initial  conditions  are  A,  =  lat  x  =  0.  Because  F  >  0,  A. will  eventually  go  to  0  within  a 
finite  time,  Tend.  At  this  terminal  time,  all  the  droplets  will  have  been  converted  to  gas.  Tend  is 
determined  by  integrating  equation  (45)  taken  between  the  limits  of  X  =  1  to  X  =  0: 


"Fend  —  J 


dV 


W,Q0) 


(47) 


Once  X(t)  is  obtained,  we  determine  p(x)  from  equation  (30)  and  (41). 


rDo  dA, Tj-jq  dA<  dx 

A  dt  A  dx  dt 


Pobf(mt)A) 


(48) 


and 


p(x)  =  y  (^go^poXl  aoL(1  +  KT(p0,T|DO,g0,l-A.3))  .  (49) 

Po  (l-y(A.3,g0>TiDoXM 


In  summary,  the  use  of  dimensional  analysis,  combined  with  specific  knowledge  of  the 
differential  equations  that  define  the  system,  has  reduced  the  solution  of  a  complicated 
burning-rate  problem  to  the  evaluation  of  a  single  integral.  Only  the  following  four 
dimensionless  parameters,  that  is, 


rbo=nvDO,  g0=  — ,  ct0  =  bpG0 , and  |30  = 


Pl 


c  T„ 
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are  required  to  distinguish  one  solution  from  another.  Originally,  21  parameters  were  used  to  set 
up  the  problem. 

3.  Lumped  Parameter  Model  for  RLPG 

3.1  Scaling  of  Differential  Equations.  The  purpose  of  this  section  is  to  use  the  direct 
assault  technique  of  dimensional  analysis  theory  in  order  to  gain  some  insight  into  how 
dimensionless  parameters  might  be  useful  for  understanding  pressure  oscillations  caused  by 
combustion  pressure  instabilities.  The  model  we  have  selected  here  is  a  simplification  of  the 
lumped  parameter  RLPG  configuration  originally  developed  by  Coffee  [3].  This  model  is 
physically  and  mathematically  not  capable  of  representing  pressure  oscillations.  Its  primary  use 
in  this  study  is  to  obtain  a  preliminary  assessment  of  the  potential  gains  in  simplification  using 
dimensional  analysis  theory.  Pressure  oscillations  are  considered  in  section  4. 

The  basic  notion  was  that,  if  dimensional  analysis  was  not  useful  in  predicting  scaling 
relationships  in  this  major  simplification  of  an  RLPG  system,  it  was  not  likely  to  work  in  a  more 
realistic  model.  Fortunately,  scaling  relationships  derived  for  this  model  were  confirmed  by 
numerical  analyses  using  the  ARL  code.  These  test  cases  are  presented  at  the  end  of  this  section. 

We  show  that  a  combination  of  analytical  techniques,  combined  with  the  direct  assault  on  the 
resulting  equations,  is  very  useful.  The  mathematical  foundation  for  the  study  of  pressure 
oscillations  makes  three  simplifications  to  Coffee’s  model:  the  piston  is  weightless,  the  liquid 
propellant  is  incompressible,  and  the  projectile  is  motionless.  The  first  two  assumptions  were 
found  quite  acceptable  using  numerical  evaluations  of  previous  studies.  The  third  assumption  is 
equivalent  to  assuming  an  infinitely  heavy  projectile. 

The  physical  model  is  similar  to  that  shown  in  Figure  2,  an  in-line  shower-head  RLPG  [3]. 
In  this  diagram,  A!  is  the  area  of  the  piston  on  the  liquid  side,  A3  is  the  area  of  the  piston  on  the 
chamber  side,  and  Av  is  the  area  of  the  vents  that  permit  propellant  to  flow  from  its  reservoir 
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Figure  2.  In-Line  Shower-Head  RLPG. 


into  the  combustion  chamber.  V3  is  the  time-dependent  volume  of  the  combustion  chamber,  and 
V,  is  the  time-dependent  volume  of  the  propellant  reservoir. 

Coffee  [3]  has  developed  a  complete  set  of  equations  that  describes  the  system  in  the  absence 
of  the  three  simplifications  made  here.  We  do  not  duplicate  his  analysis  in  detail.  We  do, 
however,  include  a  limited  number  of  equations  that  have  led  us  to  the  final  result. 

The  net  result  of  our  analysis  is  the  development  of  the  dimensionless  differential  equations, 
equations  (72)  and  (73).  These  equations  allow  us  to  make  scaling  predictions  that  include  time 
behavior  based  on  dimensional  analysis. 

As  the  pressure  in  the  chamber  increases,  the  piston  is  pushed  to  the  left.  If  Vio  is  the  initial 
volume  of  the  liquid  propellant  chamber,  the  volume,  V1?  at  time  t  is  given  by 

Vi  =  V10  -  SpsAj.  (50) 

Sps  is  the  displacement  from  the  initial  condition,  measured  positive  when  moving  to  the  left. 
It  is  related  to  the  piston  velocity,  vps ,  via  the  equation 


dS, 

dt 


=  V 


ps  ’ 


(51) 


vps  is  also  positive  when  the  motion  is  to  the  left. 

When  the  piston  is  weightless,  there  is  no  net  force  on  the  piston.  The  liquid  and  gas  forces 
on  the  piston  must  be  balanced,  giving  the  equation 


p3  (A3  -  Av)  =  pi  (Ai  -  Av). 


(52) 
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Here,  p3  is  the  gas  pressure  in  the  chamber  and  pi  is  the  liquid  pressure  in  the  reservoir.  Liquid 
enters  the  chamber  from  the  reservoir  through  the  vents.  Equation  (52)  immediately  provides  the 
relationship  between  the  differential  pressure  and  the  pressure  in  the  combustion  chamber.  We 
have 


P1-P3  =TIP3> 


(53) 


where 


3  At 

Ai  —  Av 


(54) 


Since  A3-Av  >A,-AV ,  the  pressure  in  the  reservoir  is  greater  than  that  of  the  chamber  and 
liquid  will  flow  into  the  chamber. 

The  time  dependence  of  the  mass  flow  rate,  m ,  into  the  chamber  is  controlled  only  by  the 
time  dependence  of  p3 .  This  is  easily  shown  starting  with  the  equation 

m  =  CdpLAvv3,  (55) 


where  Cdis  the  discharge  coefficient  and  v3is  the  velocity  of  the  liquid  as  it  enters  the  chamber. 
The  specific  value  of  Cd  is  inconsequential  for  our  analysis,  and,  for  this  reason,  it  is  set  equal 
to  1.0.  Equation  (55)  then  becomes 


m  =  pLAvv3, 


(56) 


where  v3  is  determined  from  Bernoulli’s  law  for  streamline  flow.  The  general  expression  is 
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(57) 


2  _  2  ,  2(Pi  p3) 

v3  —  V1 


Assuming  that  v,  =  0  and  inserting  the  resulting  expression  into  equation  (56)  gives 


m  =  Av>/2pL(p1-P3)  =Av>/2pL'np3  =WpT’  (58) 

where 


\|;  =  AvA/2pjJ.  (59) 

Using  equation  (58)  in  combination  with  the  assumption  of  incompressibility  of  the  liquid 
enables  us  to  relate  the  mass  flow  rate  to  the  velocity  of  the  piston.  The  decrease  in  volume  of 
the  reservoir  is 


d % 

dt 


=  ~VPSA1- 


(60) 


With  vpsAjPL  as  the  mass  loss  in  the  reservoir  per  unit  time,  it  must  equal  the  mass  flow  rate 
into  the  chamber.  This  gives 


VpsAl=AvV3-  (61) 

The  temporal  relationship  between  the  mass  density,  p3,  and  volume,  V3,  within  the 
chamber  is  derived  in  terms  of  the  initial  conditions  from  equations  (58)  to  (61).  The  increase  in 
chamber  volume  is  given  by 

^■=vpA,=|^7.  (62) 
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where 


g- A3Av  (63) 

Ai  V  Pl 

The  increase  in  total  mass  (equation  [58])  of  gas  in  the  chamber,  M3 ,  is  given  by 


dMj 

dt 


(64) 


From  equations  (62)  and  (64),  we  obtain 


dV3  _  a 
dM3  "  pL  ’ 


(65) 


where 


a  =  >  1 .  (66) 

A, 

Integrating  equation  (65)  between  the  initial  conditions  (denoted  by  the  additional  subscript  “0”) 
and  the  state  at  a  later  time  gives 


1- 


v3=v307 


qp3o 

Pl 


1- 


CCp, 

Pl 


Using  equation  (67)  in  conjunction  with  the  Noble-Abel  equation  of  state, 


(67) 


24 


_  (Y-l)p3CvT3 
3  (l-bp3) 


(68) 


and  the  energy  balance  equation  derived  by  Coffee  [3], 


P3C1 


dT3 

dt 


91np3 

3lnT, 


ml 


eL  + 


(l  +  r|)p3 
Pl 


CpT3  bp  3 


/P3 


(69) 


reduces  the  solution  of  the  problem  to  two  simultaneous  time-dependent  equations  for  the 
pressure  and  density.  The  details  are  tedious  but  straightforward.  From  these  two  equations,  we 
can  derive  two  equivalent  dimensionless  equations.  The  variables  in  these  dimensionless 
equations  are  the  dimensionless  density,  p ,  and  the  dimensionless  pressure,  p .  For  this  transient 
problem,  it  is  convenient  to  define 


and 


P3  = 


P3 

P30 


(70) 


(71) 


The  derived  time-dependent  equations  are 


dP3  _  (1-PP3)2 

dt  (1-3) 


(72) 


and 
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(73) 


dp 3  yp3  dp3  |  (l-PMY-l)(hi 

"^"-P3(l-kp3)  dt  (1  —  P)(l  —  kp3 ) 


Here,  t  is  a  dimensionless  time  and  is  defined  by  the  equation 


(74) 


where  x  is  a  characteristic  time  constant  for  the  system;  it  is  defined  by  the  equation 


1_  _  IPsoPl  (V^l)^-v 
t“V  P30  V30 

All  the  variables  and  parameters  in  equations  (72)  and  (73)  are  dimensionless.  They  are 


hj  =e!+gp3, 


P3(7~kP3> 

P3(7-l) 


•  g  =  —  (1+T1), 
Pl 


•  k  =  bp30,  and 
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•  B  =  ^3  ^30- ,  where 
AiPl 


2  _  P30 
P30 


C30  - 


A  close  examination  of  equation  (73)  and  the  associated  definition  shows  that  the  time 
behavior  for  the  dimensionless  pressure  and  density  as  a  function  of  dimensionless  time  are 
function  of  only  five  dimensionless  parameters.  These  parameters  form  the  set,  Q,  where 


O— /hn  A3  P30  Av  _£i_, 

i2  — \bp30»  9  9  9  2  *  * 

Ai  Pi  Aj  c30 


(76) 


Systems  that  have  the  same  values  of  the  five  dimensionless  constants  in  equation  (76)  will 
exhibit  identical  time  behavior  as  a  function  of  dimensionless  time,  t.  If  systems  have  the  same 
values  of  equation  (76)  but  differ  in  the  characteristic  time  constant,  x ,  their  motion  in  real  time 
will  be  self-similar  or,  simply,  similar;  that  is,  the  motion  will  look  identical  but  in  different  time 
scales. 

We  have  numerically  confirmed  the  self-similarity  predictions  of  equations  (72)  and  (73)  for 
a  variety  of  test  cases  using  the  ARL  computer  model.  The  success  of  this  analysis  paved  the 
way  for  the  study  of  the  more  sophisticated  Sandia  National  Laboratories  (SNL)  test  fixture 
discussed  in  section  4. 

3.2  Formal  Dimensional  Analysis  Approach.  Table  3  lists  the  variables,  parameters,  and 
dimensionality  for  the  model  described  in  section  3.1.  For  convenience  in  constructing  the 
II-functions,  we  have  divided  the  variables  in  classes  that  contain  the  same  dimensions.  It  is  of 
interest  to  see  how  the  formal  theory  of  dimensional  analysis,  based  on  Buckingham’s 
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Table  3.  Variables,  Parameters,  and  Dimensionality  for  the  Shower-Head  RLPG 
Model  Described  in  Section  3.1 


Symbol 

Definition 

Dimensions 

Index 

Class 

Time 

t 

Time 

T 

Ul 

Basis 

Areas 

A! 

Piston  Area  on  Liquid  Side 

L2 

U2 

Basis 

As 

Piston  Area  on  Combustion  Side 

L2 

U5 

1 

Av 

Vent  Area 

L2 

U6 

1 

Volumes 

v30 

Initial  Volume  of  Chamber 

L3 

U7 

2 

v3 

Instantaneous  Volume  of  Chamber 

L3 

u8 

2 

Density  Related 

Pi 

Liquid  Density 

ML’3 

U9 

3 

P30 

Initial  Gas  Density 

ML’3 

UlO 

3 

P3 

Instantaneous  Gas  Density 

ML’3 

Uil 

3 

b 

Covolume 

M^L3 

u3 

Basis 

Pressures 

Pi 

Fluid  Pressure 

ML’lrr2 

Ul2 

4 

P30 

Instantaneous  Pressure  in  Chamber 

ML’lrr2 

Ui3 

4 

P3 

Initial  Pressure  in  Chamber 

ML’lrr2 

Ul4 

4 

Velocities  | 

v3 

Propellant  Velocity  Into  Chamber 

LT1 

Ui5 

5 

VPS 

Piston  Velocity 

LT1 

Ul6 

5 

Energetics 

Cp 

Specific  Heat  at  Constant  Pressure 

L2T"2e’i 

Ul? 

6 

6c 

Chamber  Temperature 

e 

u4 

Basis 

Propellant  Energy  Per  Unit  Mass 

l2t~2 

Ul8 

7 
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II-theorem,  determines  the  dimensionless  variables.  In  the  next  section,  we  show  the 
equivalence  between  the  two  methods. 

For  a  system  involving  “m”  dimensions  and  a  total  of  “n”  variables,  the  n-functions  are 
labeled  and  constructed  according  to  the  following  scheme: 

n n_m  =  u *n-m u  2n'm u 3n~m u  u n ;  n  > m .  (77) 

By  expressing  the  un  in  terms  of  their  dimensions,  we  determine  the  IIn-m  by  requiring  the 
product  of  equation  (77)  to  be  nondimensional.  Those  un  that  have  the  same  dimensions  will 
have  the  same  values  of  wn_m,  xn_m,  yn_m,  and  zn_m. 

In  our  case,  there  are  m  =  four  dimensions,  M,  L,  T,  and  0,  and  a  total  of  n=18 
variables/parameters.  There  will  then  be  18-4=14  II-functions.  The  basis  functions  are 

•  ui  =  t, 

•  u2  =  Ai, 


•  u3  =  b,  and 

•  u4  =  0C. 

Working  through  the  procedure  developed  in  section  2.2,  we  determine  the  14  II-functions. 

•  Class  1: 


(1)  ni  =tw' Ai‘byi0*'A3  =  Tw‘ L2x‘ M~y‘ L”3y‘ 0Z' L2  =—L  and 

^■i 
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A 

(2)  n2  =  tw‘Ax‘byi0ZlAv  =  T  "2  L2*2  M ”y2 L_3y2 0 Z2  L2  =-JL,  where 

A1 

w,  =  w2  =0;  x !  =  x 2  =-l;  y!  =  y2  =°;  zi  =  z2  =0- 

•  Class  2: 

(3)  n3  =tW3AX3by30Z3V3O  =  T  W3  L2*3  M _ys L~3y3  0 23  L3  =  -^-  and 

Ai 

V 

(4)  n4  =tW4AX4by40Z4V3  =TW4L2x4M'y4L3y40ZlL3  =—575-,  where 

A1 

w3  =  w4  =0;  x3  =  x4  =  -3/2;  y3  =  y4  =0;  z3  =  z4  =0. 

Following  the  same  procedure,  we  state  the  results  for  the  remaining  IT-functions  in  simpler 
form  by  deleting  references  to  the  values  for  w,  x,  y,  and  z. 

•  Class  3 


(5)  n5  =  twAxby0zPL  =  TwL2xM'yL3y0zML-3  =  bpL , 

(6)  n6  =  tw Axby0zp3O  =  TwL2xM"yL3y0zML"3  =  bp30 ,  and 

(7)  n7  =  tw A1xby0czp3  =  TwL2xM"yL3y0zML'3  =  bp3 . 


•  Class  4: 


(8)  n8  =  tw A^by0zPi  =  TwL2xM_yL3y0zML-1T-2  =  t'A^Pi, 

(9)  n9  =  tw A^by0zp3O  =  TwL2xM”yL3y0zML”1T_2  =  t^bpjo ,  and 

(10)  ni0  =twAxby0zp3  =  TwL2xM"yL3y0zML_IT“2  =t2Aj1bp3. 
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Class  5: 


(11)  nu  =twA^by0czv3  =TwL2xM_yL3y0zLT"1  =tA'1/2v3  and 

(12)  ni2  =  twAxby0zvps  =  TwL2xM“yL3y0zLT_1  =  tA‘1/2vps . 


•  Class  6: 


(13)  n]3  =  twAxby0zCp  =TwL2xM-yL3yezL2T"2e-1  =t2A“16cCp. 


•  Class  7: 


(14)  ni4  =twAxby0zeL  =TwL2xM_yL3y0zL4T“4  =t2A^eL. 


The  general  solution  to  problem  is  given  by 


F(IIi,  n2...rii3,  II14)  =  0. 


(78) 


Because  of  the  large  number  of  IT-functions  and  their  cumbersome  forms,  equation  (78)  is  not 
especially  useful.  It  is  necessary  to  reduce  both  the  number  and  complexity  of  the  II-functions. 
We  now  discuss  methods  for  accomplishing  this. 


It  will  be  recalled  that,  in  the  direct  assault  approach  of  the  previous  section,  we  ended  up 
with  two  ordinary  time  dependent  differential  equations  for  the  dimensionless  pressure, 
p3  =  p3  /p30,  and  the  dimensionless  density,  p3  =  p3  /p30.  Without  even  solving  equations  (72) 
and  (73)  of  the  previous  section,  we  can  write  their  solutions  in  the  forms 


P3 


%(-,bp 

P  T 


30  > 


A3 

A: 


P30  Av  eL 
»  .  >  2  ' 

Pi  A,  c30 


(79) 
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and 


A3  p 


a  iTf  /  L  xjL3  K30  *  *v  \ 

P3  ”  ’^P30»  A  >  ’  A  ’2^  ’ 


Ai  P i  A!  c 


(SO) 
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with  'Pp  and  %  as  functions  of  their  respective  arguments. 

It  is  recalled  that  equations  (72)  and  (73)  are  two  coupled  equations.  It  is  possible  to  solve 
these  equations  in  quadrature  by  initially  formulating  the  single  differential  equation. 


d^ 

dp3 


:f(P3>P3>a)- 


(81) 


The  result  of  the  foregoing  discussion  is  that  the  solution  of  equations  (72)  and  (73)  can  be 
expressed  in  the  form 


rjt  *  v  A3  P30  Ay  ®L_\  _  /-) 

F'CPs’Ps’bPso’  A  »  ’  *  ’2 

A1  Pi  A1  C30 


(82) 


where  F  is  a  function  of  the  seven  variables  of  equation  (82). 

By  comparing  equations  (78)  and  (82),  we  immediately  see  that  the  former  set  contains  twice 
the  number  of  variables  than  the  latter.  This  condition  resulted  by  solving  some  of  the  algebraic 
equations  in  the  system.  That  is,  we  were  able  to  eliminate  some  of  the  variables  altogether  (e.g., 
temperature,  piston  velocity,  etc.)  by  expressing  them  in  terms  of  other  variables.  It  is  apparent 
that,  if  we  are  able  to  formulate  all  the  equations  to  begin  with,  as  we  can  in  this  simple  problem, 
there  is  probably  no  need,  or  very  limited  need,  for  using  formal  dimensional  analysis. 

3.3  Equivalence  Between  Approaches.  It  is  of  interest  to  compare  equation  (82)  with  the 
14  11-functions  previously  generated.  For  this,  we  need  to  recall  that  the  II-functions  are 
members  of  a  group,  and  satisfy  the  following  group  properties. 
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•  If  Ili  is  a  member  of  the  group,  so  is  11“ ,  where  n  is  any  positive  or  negative  real  number. 

•  If  III  and  Ilj  are  members  of  the  group  so  is  II”  11“ ,  where  n  and  m  are  any  positive  or 
negative  real  numbers. 

Thus,  there  is  nothing  unique  about  the  14  previously  generated  II-functions.  Using  the  two 
group  properties  just  listed,  a  completely  different  set  of  II-functions  can  be  constructed,  which 
is  better  suited  to  physical  interpretation.  For  example,  the  pervasive  appearance  of  time,  t,  in  so 
many  of  the  II-functions  does  not  enhance  one’s  physical  insight  into  the  nature  of  the  solution. 
The  dependence  of  so  many  of  the  II-functions  on  t  is  easily  rectified  using  the  group  properties. 

For  now,  all  we  want  to  do  is  to  show  how  equation  (82)  can  be  expressed  in  terms  of  a 
modified  subset  of  the  14  original  II-functions  using  the  group  properties.  The  first  few  steps  are 
easy  to  see: 


kP30  —  Hg, 


P30  _  Hg 

p,  n5  ’ 


(83) 

(84) 

(85) 

(36) 

(87) 


and 
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(88) 


P3  _  ^7 
P30 


Now,  consider  the  term  eL  /  c30 .  We  express  c30  as  follows: 


2  _  P30  _  ^9 

C30“  „  -.2,  -It 


P30  ^  ^-1  ^5  • 


Dividing  II14  by  (II*)2  gives 


eL  _  ^14^-6 

c2  n 

ton  IAq 


Since  there  is  nothing  unique  about  any  of  the  Ilj,  the  ratios  of  II-functions  depicted  in 
equations  (87)  to  (90)  are  just  as  valid  as  the  original  set  of  II-functions.  We  now  define  a  new 
set  of  II-functions  (using  a  “bar”  over  the  II  symbol)  in  terms  of  the  initial  set: 


n,=np  n2=n2,  n3=n3, 


n„  =  n4,  n.=5i, 


n6  ~ri6,  n7  —  ,  n8— n8, 


rin  _  Hu  *  1^12  1^12 »  1^13  H13 » 


tt  -11  n  = 

Il9—  il9,  iA10  J-J.  » 


TT  —  ILiHl 

14  ■  nD 


Equation  (82)  now  becomes 


p  (nltn2fn5,n6,n7,n10,n14)  —  o . 
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The  foregoing  equation  or  its  equivalent,  equation  (82),  is  actually  easier  to  use  than  it  would 
appear.  The  reason  for  this  is  that  some  of  then -functions  are  actually  constants,  or  parameters 
that  do  not  change  much.  For  example,  II5  and  II6are  expected  to  be  nearly  constant  under  most 

conditions.  If  the  combustion  chamber  is  designed  to  retain  constant  area  ratios,  then  III  and 
n2  are  constants  and  equation  (82)  becomes 


F,(P3^3>%=0-  (93) 

C30 

Equation  (93)  shows  that,  under  the  aforementioned  conditions,  the  relationship  between 
pressure  and  density  depends  only  on  the  single  dimensionless  variable,  eL  lc\0 .  The  time  scale 
associated  with  the  behavior  of  pressure  and  density  is  determined  from  equations  (72)  and  (73). 

For  systems  where  there  is  high  confidence  that  the  physics  have  been  accurately  modeled, 
the  direct  assault  approach  is  probably  the  best  method  for  determining  the  significant 
dimensionless  variables.  On  the  other  hand,  when  the  equations  that  describe  the  system  are  not 
well  defined,  the  formal  theory  can  provide  guidance  in  deducing  significant  dimensionless 
parameters.  The  virtue  of  being  able  to  solve  a  portion  of  the  problem  is  the  reduction  of  the 
number  of  II-functions. 

4.  Dimensional  Analysis  for  a  Test  Fixture 

4.1  Physical  Model.  A  simplified  problem  has  been  chosen.  A  test  fixture,  Figure  3,  was 
developed  at  SNL  to  study  pressure  oscillations  [15].  The  combustion  chamber  in  the  fixture  is 
basically  a  cylinder  with  a  1-in  radius.  Liquid  propellant  is  injected  through  a  circular  opening 
on  the  center  line  at  the  left  of  the  chamber  by  a  differential  piston.  Hot  gas,  and  possibly  some 
liquid,  flow  through  a  converging  orifice  at  the  right  into  atmospheric  pressure.  Pressure 
oscillations  are  quickly  established  in  the  chamber.  Normally,  there  is  a  single,  clear  frequency 
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corresponding  to  the  first  radial  mode  of  the  cylinder.  This  system  is  very  similar  to  a  liquid 
propellant  rocket. 

The  numerical  model  is  simplified  as  much  as  possible  for  this  exercise.  The  chamber  is 
assumed  to  be  exactly  a  cylinder,  ignoring  the  shape  of  the  converging  nozzle.  The  liquid  is 
injected  at  a  steady  rate.  The  liquid  instantaneously  breaks  up  into  droplets  of  a  specified  size 
and  ignites.  Combustion  occurs  according  to  a  pressure-dependent  bum  rate,  ApB .  The  liquid 
and  gas  have  the  same  velocity  and  the  same  pressure  at  any  point.  Viscosity  and  turbulence  are 
ignored.  Heat  transfer  to  the  liquid  is  ignored.  The  liquid  is  assumed  to  be  incompressible 
(treated  as  compressible  in  the  code).  The  outflow  is  computed  assuming  choked  flow.  At  all 
points  within  the  fixture,  the  gas  and  liquid  are  assumed  to  move  with  the  same  velocity.  The 
numerical  simulations  quickly  reach  pseudosteady  state.  The  mean  pressure  and  outflow  rate  are 
constant,  but  there  are  oscillations  about  the  mean.  For  the  dimensional  analysis,  the  outflow 
rate  is  considered  constant  and  equal  to  the  inflow  rate. 

4.2  Eulerian  Equations.  In  this  section,  the  basic  continuity,  momentum,  and  energy 
equations  are  derived  for  the  Sandia  fixture.  For  mathematical  simplicity,  the  equations  are 
rendered  in  Cartesian  coordinates.  In  actuality,  the  system  is  described  best  in  cylindrical 
coordinates.  The  scaling  relationships  and  insights  are  the  same  for  both  coordinate  systems. 

There  are  a  large  number  of  variables  and  parameters  in  the  mathematical  model  for  the 
Sandia  fixture.  All  of  these  variables  and  parameters  are  defined  in  the  List  of  Symbols.  To 
facilitate  the  presentation,  some  of  them  are  also  defined  in  this  section  and  in  section  4.3. 

The  continuity  equation  for  the  gas  is 


3(P)  ,  3(pVx)  .  9(pVv) 

IT —y 


=  4jtpLfDron . 


(94) 
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The  right-hand  side  of  equation  (94)  is  the  volumetric  increase  in  gas  mass  density  by  burning  of 
the  liquid  drops. 

Because  the  gas  and  droplets  are  assumed  to  move  with  the  same  velocity,  we  can  use  a 
volume  average  mass  density  for  the  gas  plus  droplets.  This  is  defined  by 

Pt=P+Pd>  (95) 

where  p  is  the  mass  density  of  the  gas  (gas  mass  divided  by  volume)  and  pD  is  the  mass  density 
of  the  droplets  (liquid  mass  divided  by  volume). 

The  two  momentum  equations  are 


d(pTVx)  ,  3(PtVx)  ,  a(pTVYVX)  dp 

dt  dx  dy  dx 


and 


3(PtVy)  ,  d(pTVxVv)  !  3(PtVv)  _ 
d  t  d  x  d  y  3y 


The  energy  equation  for  the  gas  is 


d  ( 1  ,72 

_  -pV  +P„ 


\  f 
+  V* 


pV(IV2+h) 


—  eLS 


D-»G  ’ 


where 


u  =  CmT  =  internal  energy  per  unit  mass, 


(96) 


(97) 


(98) 


(99) 
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(100) 


h  =  u  +  — 
P 


= enthalpy  per  unit  mass, 


and 


V2 


=  V2+V2. 


(101) 


The  pressure  is  related  to  the  density  and  temperature  via  the  equation  of  state, 


P  = 


pgRT 

l-bpG’ 


where  pG  is  the  true  gas  density.  It  is  given  by 


Pg  — 


(102a) 


(102b) 


where  n  is  the  number  density  of  the  liquid  drops. 

We  now  consider  the  droplet  dynamics.  Because  equations  (95)  to  (98)  describe  the 
gas-plus-droplet  system,  there  is  no  additional  momentum  equation  for  the  droplets.  We  assume 
that  heat  transfer  between  droplets  and  gas  is  neglected  and  that  droplets  supply  energy  to  the  gas 
when  they  convert  to  gas.  These  assumptions  preclude  the  need  for  an  energy  balance  for  the 
droplets.  The  only  relevant  equation  for  the  droplets  is  the  mass  balance  equation.  This  equation 
is  given  by 


3(Pd)  |  ^(PdVx)  ,  a(PpVy) 
3t  dx  dy 


=  -Sd->g  =-4ttpLfDr2n, 


(103) 
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where 


4k  3 

Pd  =  ^Pl^d  —  npL  ^ 


(104) 


and 


rD  =  ApE 


(105) 


Dividing  both  sides  of  equation  (103)  by  pL  and  using  equation  (105)  gives 


3 (nr*)  t  3(nrpVx)  {  8 (nr pVY)  _  3ApBr2n  (106) 

3t  3x  3y 

Equation  (106)  is  a  partial  differential  equation  in  two  dependent  variables:  n  and  rD.  This, 
however,  can  be  reduced  to  a  single  partial  differential  equation  in  the  single  variable,  ro,  by 
recognizing  that  the  continuity  equation  for  the  number  density  is 


9(n)  4  a(nVx)  4  d (nVY )  _ Q 
3t  3x  By 


(107) 


That  is  to  say,  even  though  the  mass  of  a  droplet  is  always  decreasing,  the  droplet  always 
exists— even  though,  in  the  limiting  case,  it  may  have  zero  mass.  Differentiating  the  left-hand 
side  of  equation  (106)  and  using  equation  (107)  gives  the  following  Eulerian  description  of  the 
droplet  radius: 


2M  +  Vx  ^1+ V  =  -ApB .  (108) 

3t  3x  3y 


Equation  (108)  applies  only  for  rD  >  0. 
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In  summary,  we  have  a  set  of  eight  unknowns,  p,  pT,  Vx,  Vy,  p,  T,  n,  and  rD,  in  eight 
equations,  (94)-(98),  (102),  (107),  and  (108).  The  object  now  is  to  reduce  this  system  of 
equations  to  an  equivalent  set  of  dimensionless  equations  that  can  be  analyzed  using  a  relatively 
small  number  of  parameters. 

4.3  Self-Similar  Eulerian  Equations.  In  this  section,  we  convert  the  equations  of  the 
previous  section  to  dimensionless  forms.  We  define  the  following  set  of  dimensionless 
variables: 


p  =  p/p0,p  =  p/p0,  T  =  T/T0,n=n/n0,  rD  =rD/rD0,  Vx  =  VX/V0,and  Vy  =Vy/V0, 

where  the  baseline  values  are  denoted  by  the  subscript  0.  There  is  no  restriction  on  the  baseline 
values.  When  the  foregoing  dimensionless  variables  are  inserted  into  the  system  of  equations  of 
section  4.2,  we  obtain  a  set  of  dimensionless  equations.  The  dimensionless  continuity  equation 
is 


dp  .  p  9(pvx)  ,  p  3(pvy) 


E0npBr£, 


(109) 


where 


(110) 


p  _  VqTqo 

X  o  — 


V()rD0 

LAp0B 


(111) 


E0 


4tcpLn0rD0 

Po 


(112) 
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(113) 


Too  — 


*■00 

Apo 


and 


X  = 


(114) 


The  dimensionless  momentum  equations  are 


g(PTVx)H  p  3(&V'2 


3t 


3x 


v) ,  r  3(PtVvVx) 


8y 


—  — P  K  ^P 
"  0  °dx 


(115) 


and 


a(pTVy)  |  p  ^(pTVxVy)  ■  p 
ax  0  ax  c 


5(Ptvy)_  pr  ^P 

~W~~  0  0 ay’ 


(116) 


where 


K0 


Po 

P0V02 


and 


E 

Pt=(p+y5Fd)- 


(117) 


(118) 


The  dimensionless  energy  equation  becomes 
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dz 


Ipv’  +  PC»TT° 

2 


V2 

vo 


+P0v. 


(  A  1  _  c  TT 

pV(-V2+  m  ,  q) 
2  V2 


■\ 


+  K0P0V.(pV)  =  H0npBr2,  (119) 


where 


V 


r  d  t  3 

dx  3y 


(120) 


H0 


(121) 


and 


pT  rppT0R^ 

l-nr2(^)-bp0pl  Po  ) 

3pL 


The  dimensionless  droplet  equations  are 


(122) 


3(5)  ,  r  9(5Vx)  ,  ?  3(5Vy)  0 

3x  0  3x  0  ay 


(123) 


and 


^  +  PoVx^+PoVv^  =  -p”.  (124) 

at  ax  3y 

Simplification  of  the  dimensionless  equations  occurs  by  making  the  appropriate  selections 
for  the  baseline  parameters.  This  is  where  physical  insight  into  the  problem  is  of  great  value. 
There  is  a  temptation  to  select  values  of  pressure,  temperature,  and  density  at  their  initial  values. 
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However,  we  know  from  experience  that  within  a  broad  range  of  initial  conditions,  the  response 
of  liquid  rocket  systems  and  RLPG  designs  are  essentially  independent  of  the  initial  values  of 
pressure,  temperature,  and  density.  Under  these  conditions,  we  are  liberty  to  select 
p0 ,  T0 ,  and  p0  to  simplify  the  dimensionless  equations. 

We  select  p0  and  p0  to  make  P0  =  1  and  E0  =  1 .  This  gives 


Po  = 


fVr  V/B 

I  vOiDO 

LA 


Po  —  ^Pl^O^DO  » 


Po  (VorDO)1/B 

p0V02  3(LA)1/BpLnouD0 


(V0rDo)1/B 
3(LA)1/BpLt)D0  ’ 


and 


4 n  3 

^Ido  =  ^o^do  —  n0  2  ' 


(125) 


(126) 


(127) 


(128) 


When  the  foregoing  definitions  are  used,  the  dimensionless  time  now  becomes 


t 


TDOis  now  given  by 


(129) 


(130) 
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The  baseline  temperature,  T0 ,  is  selected  to  satisfy  the  condition 


C  T 

10  _  i 

o  A  • 

V 

vo 


This  gives 


T0  = 


(131) 


(132) 


When  the  foregoing  definitions  are  used,  the  dimensionless  time  now  becomes 


T  = 


(133) 


where 


T  —  rDQ  -  h. 
D°  Apg  V0 


(134) 


The  new  system  of  dimensionless  equations  now  becomes 


ap  a(pvx)  a(pvY)  __B_2 

_J1+.  V;_Y  =npBr^, 

Bx  By 


(135) 


a(pTVx) ,  ^(PTVx)  ,  3(PtVyVx)  k  dp 

Bx  Bx  By  °Bx’ 
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(137) 


3(pTVY)  .  3(PTVXVY)  3(ptVy)  _  „  3p 

—sT~+  3i  +_ 35  °3y 


g 

Pr=(P+Ynf^), 


(138) 


Afp(I +  f)  V V  •  fp V(- V2  +  T)  |+ K0 V  •  (pY  )  =  H0n pBfD2 , 
dx[  2  J  l  2 


(139) 


pT  r  PoT0R 

l-nr2(^-)-bp0pl  Po 

3pL 


(140) 


d(n)  ,  a(nVx)  ,  3(nVY)_0 
dx  dx  dy 


(141) 


^Oij)  ,  y  ^(ro)  |  y  d  (rp)  _  — B 

TT  s_3r  v  3y  p  ' 


(142) 


By  writing  equation  (140)  in  the  form 


1 —  ^rDTJD0  —  N0  P 


(143) 


where 


N0  =  bp0  =3bpLT|D0 


(144) 
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and 


Q0  — 


Po^pR 

Po 


W  '  *  V'Y3pL11D0VosR^ 


LA 

y  ^0rD0 


(145) 


we  arrive  at  a  system  of  eight  equations,  (135)  to  (142);  eight  dimensionless  variables, 
p ,  pT ,  Vx ,  Vy ,  p,  T,  n ,  and  rD ;  and  five  dimensionless  constant  parameters,  Ko,  Ho,  t|do,  No,  and 

Qo.  Examination  of  r|D0  and  N0  shows  that  we  can  eliminate  one  of  these  two  since  the 
covolume,  b,  is  a  constant  and  the  liquid  density  is  nearly  constant.  We  keep  r|D0  and  eliminate 
N0. 

The  solution  for  each  of  the  variables  (we  use  pressure  as  an  example)  can  be  expressed  in 
the  form 


P  —  P0fp(L  >  ^  >"LKq,H0,T|Dq,£2q)  , 


(146) 


where 


Po  = 


V  r  v 

V0rD0 

LA 


l/B 


(147) 


and  fp  is  a  function  determined  from  the  solution  of  equations  (135)  to  (142).  A  different 
f  function  applies  to  each  of  the  eight  variables. 

4.4  Confirmation  of  Scaling  Relationships.  Systems  that  have  the  same  values  of  the 
dimensionless  constant  parameters  Ko,  Ho,  T)o,  and  Qo  will  have  similar  time  and  space  behavior 
when  expressed  in  dimensionless  variables.  Consistency  between  the  current  Eulerian 
formulation  and  the  Lagrangian  formulation  used  in  the  code  allows  the  examination  of  pressure 
oscillations  to  be  made  in  the  relatively  limited  dimensional  space  of  only  four  parameters. 
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While  the  dimension  space  of  the  problem  has  been  reduced,  the  relationships  between  the 
variables  are  still  extremely  complicated.  However,  analysis  of  the  results  indicates  an 
interesting  scaling  relation.  Suppose  the  dimensions  of  the  chamber  (radius,  length,  radius  of 
inlet,  radius  of  outlet)  are  increased  by  a  constant  factor.  If  either  the  radius  of  the  injected 
droplets  is  increased  by  the  factor  or  the  pre-exponential  term,  A,  is  decreased  by  the  factor,  the 
normalizing  factor,  po,  is  unchanged.  The  nondimensional  tune,  x,  includes  the  constant,  L,  so 
the  time  scale  will  decrease  by  this  factor.  The  other  nondimensional  quantities  in  the  function, 
fp  ,  are  unchanged.  The  injection  velocity  is  left  unchanged,  so  the  mass  flow  increase  is 
proportional  to  the  vent  area  increase.  The  pressure  should  be  the  same,  except  for  the  change  in 
the  time  scale. 

This  scaling  law  for  the  pressure  seems  reasonable  physically.  If  the  dimensions  are 
increased,  pressure  waves  will  take  longer  to  reach  the  walls  of  the  chamber  and  return  to  the 
combustion  site  near  the  injector.  If  the  combustion  rate  is  decreased  proportionally,  the  pressure 
waves  will  receive  the  same  amount  of  energy  at  each  iteration. 

To  check  this  result,  numerical  simulations  were  performed  using  the  2-D  RLPG  code.  The 
baseline  simulation  for  this  fixture  assumes  the  jet  breaks  up  into  200-pm-diameter  droplets.  A 
pressure  trace  from  the  simulation  at  the  top  wall  is  shown  in  Figure  4.  Very  large  oscillations 
are  generated.  The  Fourier  transform  is  shown  in  Figure  5.  The  largest  frequency  is  a  first  radial 
mode  (around  26  kHz).  With  the  chamber  represented  as  a  cylinder  (ignoring  the  actual  nozzle 
shape  at  the  exit)  there  is  a  first  longitudinal  mode  (around  10  kHz).  There  are  also  smaller 
overtones. 

Simulation  was  performed  again  with  all  the  physical  dimensions  doubled.  The  size  of  the 
injected  droplets  is  also  doubled.  The  pressure  at  the  middle  of  the  top  wall  is  shown  in  Figure  6. 
Note  that  the  time  scale  has  been  changed.  The  Fourier  transform  is  shown  in  Figure  5.  The 
magnitudes  are  almost  the  same  as  the  baseline  case,  while  the  frequencies  are  cut  in  half.  In 
Figure  7,  the  frequencies  of  the  doubled  case  are  cut  in  half.  There  is  no  noticeable  difference 
from  the  baseline  case  when  the  time  scaled  is  thus  adjusted.  The  same  answer  is  obtained  if  the 
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Figure  4.  Sandia  Test  Fixture  -  Baseline:  Simulation  With  200-)jm-Diameter  Droplets. 
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Figure  5.  Baseline  FFT  (Line)  and  All  Physical  Dimensions  Doubled  FFT  (Dot). 
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Figure  6.  Sandia 
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Times  2  (Dot). 
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pre-exponential  is  cut  in  half  and  the  droplet  diameter  is  left  at  200  pm.  The  model  produces 
almost  exactly  the  answer  predicted  by  the  dimensional  analysis. 

This  result  only  holds  if  all  the  conditions  are  applied.  For  instance,  suppose  all  the 
dimensions  are  doubled  except  for  the  length,  which  is  left  the  same.  The  longitudinal  mode,  as 
expected,  has  a  different  frequency.  But  the  radial  mode  is  also  very  different  in  amplitude. 
Also,  if  the  dimensions  are  doubled  but  the  droplet  diameter  is  left  at  200  pm,  a  completely 
different  answer  is  computed. 

Standard  correlations  for  injected  droplet  size  show  that  the  size  depends  primarily  on  the 
injection  velocity  and  the  ratio  of  the  gas-to-liquid  density  [5].  So,  if  a  scale  up  of  the  test  fixture 
were  constructed,  we  would  expect  that  the  injected  droplet  diameter  would  be  about  the  same. 
This  implies  that,  to  actually  construct  a  scale  of  a  fixture  that  has  the  same  magnitude  of 
pressure  oscillations  as  the  original  fixture,  either  the  droplet  diameter  must  be  controlled  or  the 
burning  rate  of  the  propellant  must  be  modified. 

5.  RLPG  Considerations 

The  dimensionless  equations  for  a  gun  simulation  have  not  yet  been  derived.  The  problem  is 
essentially  transient  and  much  more  complicated  than  the  aforementioned  test  fixture.  However, 
based  on  physical  intuition,  some  simulations  were  made  to  check  whether  the  scaling  mle  still 
holds  for  gun  conditions. 

Figure  8  shows  the  pressure  at  the  middle  of  the  top  wall  of  the  chamber  for  a  30-mm  gun 
simulation.  The  injected  droplet  diameter  is  now  a  function  of  time  rather  than  a  constant.  The 
dimensions  were  then  doubled  (60-mm  gun),  and  the  injected  droplet  diameter  profile  was  also 
doubled.  Figure  9  shows  the  result.  Again,  note  the  different  time  scale.  The  curves  are  very 
similar.  The  Fourier  transforms  are  almost  the  same  for  any  time  window  during  the  firing  cycle. 
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Figure  9.  6' 


6.  Intact-Core  Model 


So  far,  all  the  simulations  shown  have  assumed  that  the  injected  liquid  instantaneously  breaks 
up  into  droplets.  This  is  a  simplification  of  a  much  more  complicated  process.  A  number  of 
other  jet  breakup  models  have  also  been  considered  [16].  The  most  useful  has  been  to  assume 
that  most  of  the  liquid  is  in  a  noncombusting  intact  core.  Droplets  stripped  off  the  core  are  small, 
on  the  order  of  10  pm.  The  length  of  the  intact  core  depends  on  the  ratio  of  the  gas-to-liquid 
density  and  the  diameter  of  the  injector.  So,  the  intact-core  length  should  scale  with  the  size  of 
the  fixture. 

A  number  of  simulations  were  performed  using  the  intact-core  model.  Consider  a  30-mm 
gun  simulation.  The  length  of  the  intact  core  and  the  size  of  the  small  combusting  droplets 
outside  the  core  vary  through  the  firing  cycle.  This  case  was  then  scaled  up  to  a  60-mm  gun.  If 
the  length  of  the  intact  core  and  the  size  of  the  drops  outside  the  core  are  doubled,  the  scaling  law 
still  holds.  The  pressure  oscillations  are  the  same  magnitude,  with  the  time  scale  doubled. 
However,  if  the  small  drops  are  kept  the  same  size  as  in  the  30-mm  simulation,  the  solution  is 
very  different  in  character.  It  was  expected  that,  once  the  droplets  became  small,  the  exact  size 
would  not  be  important.  This  turned  out  not  to  be  the  case.  Thus,  the  scaling  law  for  the 
intact-core  model  is  essentially  the  same  as  for  the  droplet  model. 


7.  Conclusion 


Dimensional  analysis  has  been  applied  to  the  equations  governing  RLPGs.  A  scaling  law 
that  relates  fixtures  of  different  sizes  has  been  derived.  The  numerical  model  agrees  almost 
exactly  with  the  scaling  law.  This  increases  our  confidence  both  in  the  dimensional  analysis  and 
the  numerical  accuracy  of  the  computational  fluid  dynamics  code. 

It  is  useful  to  be  able  to  do  experiments  in  a  smaller  scale  fixture.  Unfortunately,  the  scaling 
rule  derived  here  indicates  that  either  the  injected  droplet  size  or  the  bum  rate  also  has  to  be 
modified.  The  droplet  size  is  not  under  our  control  for  practical  gun  injectors,  and  the  droplet 
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size  is  expected  to  be  about  the  same  for  fixtures  that  operate  at  similar  pressures.  There  has 
been  some  work  in  developing  propellant  with  a  faster  bum  rate  [17].  However,  the  bum  rate 
cannot  be  controlled  with  any  precision. 

It  has  long  been  noticed  that  pressure  oscillations  in  155-mm  guns  are  very  different  in 
character  than  oscillations  in  30-mm  guns.  Even  when  the  same  modes  occur,  the  magnitudes 
are  different.  The  present  work  indicates  that  this  should,  in  fact,  be  expected.  While  we  can 
scale  the  physical  dimensions  of  the  30-mm  gun,  we  cannot  control  the  droplet  size  or  the  bum 
rate  of  the  propellant. 
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